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Abstract 



Recently one of us derived the action of modified gravity consistent with the 
holographic and new-agegraphic dark energy. In this paper, we investigate the sta- 
bility of the Lagrangians of the modified gravity as discussed in [M. R. Setare, Int. 
J. Mod. Phys. D 17 (2008) 2219; M. R. Setare, Astrophys. Space Sci. 326 (2010) 
27]. We also calculate the statefinder parameters which classify our dark energy 
model. 
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1 Introduction 



Nowadays it is strongly believed that the universe is experiencing an accelerated expan- 
sion. Recent observations from type la supernovae [Ij in associated with Large Scale 
Structure [2J and Cosmic Microwave Background anisotropics |3] have provided main evi- 
dence for this cosmic acceleration. There are two ways to explain the current accelerated 
expansion of the universe. The first one is to introduce some unknown matter, which is 
called dark energy in the framework of general relativity. 

On the other hand the nature of dark energy is ambiguous. The simplest candidate 
of dark energy is a cosmological constant with the equation of state parameter u = —1. 
However, this scenario suffers from serious problems like a huge fine tuning and the 
coincidence problem [3]. Alternative models of dark energy suggest a dynamical form 
of dark energy, which is often realized by one or two scalar fields. In this respect, dark 
energy components such as quintessence |5], K-essence [6], tachyon [7j, phantom [8], ghost 
condensate and quintom [9J, and so forth. 

Although going beyond the above effective description requires a deeper understanding 
of the underlying theory of quantum gravity (TU] unknown at present, physicists can 
still make some attempts to probe the nature of dark energy according to some basic 
quantum gravitational principles. An example of such a paradigm is the holographic 
dark energy scenario, constructed in the hght of the holographic principle [11], [121 [131 E] • 
Its framework is the black hole thermodynamics [20] and the connection (known from 
AdS/CFT correspondence) of the UV cut-of of a quantum field theory, which gives rise 
to the vacuum energy, with the largest distance of the theory [TT]. Thus, determining an 
appropriate quantity L to serve as an IR cut-off, imposing the constraint that the total 
vacuum energy in the corresponding maximum volume must not be greater than the mass 
of a black hole of the same size, and saturating the inequality, one identifies the acquired 
vacuum energy as holographic dark energy: 

with G the Newton's gravitational constant and c a constant. The holographic dark 
energy scenario has been tested and constrained by various astronomical observations 
[T5| [T6l [T71 [T8l [T9] and it has been extended to various frameworks [211 l22l |23] . 

More recently a new dark energy model, dubbed agegraphic dark energy has been 
proposed. These models take into account the Heisenberg uncertainty relation of quantum 
mechanics together with the gravitational effect in general relativity. The agegraphic dark 
energy models assume that the observed dark energy comes from the spacetime and matter 
field fluctuations in the universe [231 [221 [26]. Since in agegraphic dark energy model the 
age of the universe is chosen as the length measure, instead of the horizon distance, the 
causality problem in the holographic dark energy is avoided. The agegraphic models of 
dark energy have been examined and constrained by various astronomical observations 

Another alternative approach to explain the universe's late-time acceleration is mod- 
ifying the General Relativity itself [2H], and in the simplest case replace R with f{R) 
in the action which is well known as f{R) gravity. Here f{R) is an arbitrary function 
of scalar curvature. Although there are some works with related subjects on crossing of 
the phantom divide line in the framework of modified gravity [211 [30] , but Ref . [31] was 



2 



the first paper that has investigated a modified gravity model realizing oj across -1. The 
authors of Ref. |3T] have shown an explicit model of modified gravity in which a crossing 
of the phantom divide can occur and relation between scalar field theories with property 
of u crossing -1 and the corresponding modified gravity theories have been investigated. 
Observationally the possibility of phantom crossing was suggested by early supernova 
data sets [32] while the situation concerning present data is rather ambiguous [33] . 

In paper [31], using the holographic model of dark energy in spatially fiat universe, one 
author of this paper has obtained equation of state for holographic dark energy density 
in framework of modified gravity for a universe enveloped by Rh as the system's IR cut- 
off. Also he has developed a reconstruction scheme for the modified gravity with f{R) 
action. He could to obtain a differential equation for f{R), the solution of this differential 
equation give us a modified gravity action which is consistent with holographic dark energy 
scenario. In j3Sj this investigation extended to new agegraphic dark energy model. 

2 Holographic dark energy 

The energy density of holographic dark energy is given by (SttG = 1) 



Pd = Sc^Rh^, 



(2) 



where c is a constant and Rh is the future event horizon specified by 




(3) 



a{t) 



The critical energy density pcr is given by the following expression 

Per = 3if^, 



(4) 



which helps us to write the dimensionless density parameter for dark energy 




(5) 



Differentiating ^ and using ([s]), we are able to write 



Rh — HRfi — 1 



c 



(6) 




- 1 



The energy conservation equation for (any form of) dark energy is 




(7) 



Differentiating ^ w.r.t t and using (|5| yields 





Using (g in yields 




(9) 



the equation of state parameter for the holographic dark energy. 
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3 Stability of modified gravity with HDE 



It is shown in [3l] that a modified gravity consistent with the holographic dark energy in 
fiat space has the following form: 

fiR) = C,R^ + C^W + ""'^^"2°^'^ , (10) 

where is a constant which is related to the scale factor as a{t) = ao(ts —t)^°, tg and Qq 
are also constants. The expressions of u, v, Ci and C2 are discussed in [M] . 

For a phenomenological f{R) model to be viable, it must meet a certain list of viability 
conditions [36]: classical and quantum stability which involves f'{R) > and f"{R) > 
(the first positive derivative means that gravity is attractive and the graviton is not a ghost 
while the second positive derivative is used to avoid the Doglov-Kawasaki instability); 
stable Newtonian limit for all values of R; absence of deviations from general relativity 
and existence of a future de Sitter asymptote. We are interested here to check the stability 
of our f{R) models only. The stability can be checked by calculating the double derivative 



of/(i?) in (10): 



fiR) = Ciu{u - l)R--' + C,v{v - l)R''-\ (11) 

Now f"{R) > if Ci > 0, C2 > 0, M > 1 and f > 1. This case corresponds to stability 
while a negative double derivative results in instability. We are not interested in the later 
case. There are some suggestions that the instabilities associated with f{R) action can 
be removed by adding quadratic and higher order terms along with the Ricci scalar R 
[37] . Thus due to abundance of free parameters, the instabilities associated with f{R) 
action are naturally alleviated. Note that when f"{R) = for a suitable choice of model 
parameters (e.g. Ci = = C2) then the theory reduces to general relativity and is of no 
interest here. 

4 New agegraphic dark energy 

The energy density of the new agegraphic dark energy is given by 

PD = Sn^r]-\ (12) 
where n is a constant and 77 is the conformal time defined by 





Using (12) and we obtain 



Differentiating (12), we get 



Using (15) in ([T]), we obtain 



n 

HV=^- (14) 

V 



2 

pD = pD- (15) 

arj 



2 I — 

l^NADE = -l+7^ y^D- (16) 

San ^ 
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5 Stability of modified gravity with NADE 



It is shown in ^35j that a modified gravity consistent with the agegraphic dark energy in 
fiat space has the following form: 

f{R) = C,R^ + C,R^ + ^^. (17) 

The expressions of C4 and C5 are discussed in [35]. 



Now stabihty can be checked by calculating the double derivative of f{R) in (17): 



dR^ ^ ' ^ ' ^ ' ' ' ^ ^i?i+^o- 

Now f"{R) > if C3 > 0, C4 > 0, C5 > 0, M > 1, f > 1 and ho > 1. This case corresponds 
to stability while a negative double derivative results in instability. We are not interested 
in the later case. It should be noted that the order of these coefficients should of the 
order of Hq ~ 10~^^eV^, so that the correction terms play their role at the present 
time to produce cosmic acceleration. We would also mention that at present there is 
no indication what should be the correct modification of the Einstein's general relativity 
from the observations. On the other hand, some recent astrophysical observations of 
weak lensing and galaxy velocities indicate that general relativity is still the best fit of 
the astrophysical data [ID] . The present considered models contain numerous parameters 
and observational data can constrain them only. 



6 Statefinder diagnostic 



In this section, we calculate the statefinder parameters for the above two models of dark 
energy. Sahni et al [H] introduced a pair of cosmological diagnostic pair {r, s} which 
they termed as Statefinder. The two parameters are dimensionless and are geometrical 
since they are derived from the cosmic scale factor alone, though one can rewrite them 
in terms of the parameters of dark energy. Additionally, the pair gives information about 
dark energy in a model independent way i.e. it categorizes dark energy in the context 
of background geometry only which is not dependent on the theory of gravity. Hence 
geometrical variables are universal. Also this pair generalizes the well-known geometri- 
cal parameters like the Hubble parameter and the deceleration parameter. This pair is 
algebraically related to the equation of state of dark energy and its first time derivative. 
The statefinder parameters {r, s} are defined as [HI H2] 



r — 1 



(19) 



3(g-l/2)' 

Note that in the derivation of the above parameters, a spatially fiat FRW spacetime is 



assumed. A useful alternative form of (19) is 

3 



1 - -^D 

2 



H 



-1 
3ujd 



UJd_ 



- 3wd(1 + wd) 



(20) 
(21) 
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One can immediately see that for cosmological constant with constant equation of state 
{un = —1), we have {1,0}. Moreover {1,1} represents the standard cold dark matter 
model containing no radiation while Einstein static universe corresponds to {oo, —oo} [33]. 
In literature, the diagnostic pair is analyzed for various dark energy candidates including 
holographic dark energy [H] , agegraphic dark energy [35] , quintessence [16] , dilaton dark 
energy [17], Yang-Mills dark energy [18], viscous dark energy [19], interacting dark energy 
[SO] , tachyon [HI], modified Chaplygin gas [5^, f{R) gravity [S3] and dark energy model 
with variable constants [3^ to name a few. 

The dimensionless dark energy density parameter is [3H [35] 



D 



3H^ 



f{R)-Q{H + H'-H^^)f{R) 



(22) 



where prime denotes differentiation w.r.t R. The above expression can be written as 

1 



n 



D 



3H^ 



f{R) - 6{H + H')F{R) + QHFrR 



where we have used the following relations 



F{R) = f'{R), -f\R) = f^RRR = F^RR. 



(23) 



(24) 



We proceed to calculate statefinder parameters for the NADE. Substituting (23) in (16) 

2 



^NADE 



■1 + 



3\/3anH 



f{R) - 6{H + H^)F{R) + 6HF,rR. 



(25) 



We differentiate (16) w.r.t t and obtain 

2 

^NADE 



3na 



n 



D 



1 . 



(26) 



Differentiating (|23|) leads to 



-2H 



3H^ 
1 

'^3IP 



f{R) - 6{H + H')F{R) + GHF^rR 



[FR - 6iH +2HH)F + GHF^rR + QH{F^rrR? + F^rR)]. (27) 



Similarly, we obtain parameters for the HDE: 

-1 2 



UJHDE = ^ - ^y|^V/(^) - + H^)F{R) + QHF^rR. 



-1 VL 



^HDE 



D 



3c y/Tt^ 



Using (16) and (26) in (20) and (21), we obtain 

3 



rNADE 



an 



3/2 
D 



D 



an 



D 



(28) 
(29) 

(30) 



6 



snade 



3an 



_|_ ^'-D 



-3 + 



D 



Similarly we can obtain {r, s} parameters for the HDE as 



thde 



1 - -^D 
2 



Shoe 



3(i + ^x/n 



D 



+ 2(l--fiz. 



(31) 



(32) 
(33) 



7 Conclusion 

Besides compatibility with experimental data, the minimal criteria that a modified gravity 
theory must satisfy in order to be viable are [38] : reproducing the desired dynamics of the 
universe including an infiationary era, followed by a radiation era and a matter era and, 
finally, by the present acceleration epoch. Moreover, the theory must have Newtonian 
and post-Newtonian limits compatible with the available solar system observational data. 
And the theory must be stable at the classical and quantum level. 

In this paper, we have studied the stability issues associated with the HDE and NADE 
in the framework of modified gravity. Our study provides a more general framework to 
study such models. It should be noted that the models of [39j is a special case of our 
paper for the particular choices of parameters such as C3 = (or C4 = 0, but not both 
simultaneously) and u = l{v = 0) or v = l{u = 0) and C5 = —fi'^. Also note that the 
suggestion of ^9J is ruled out since it was incompatible with the observational data of 
solar system tests. A more appropriate f{R) action is provided by Setare [SU [35] and 
here we obtained constraints on the parameters of our model if the f{R) gravity theory 
is stable i.e. f"{R) > 0. The opposite case is not discussed for not general interest and 
it corresponds to Dolgov-Kawasaki instability [SB]. Finally we calculated the statefinder 
diagnostic parameters to characterize the present models of dark energy and plotted their 
trajectories in Fig 1 and 2 for a suitable choice of model parameters. The dots in both 
figures represent the LCDM model for which {r,s} = {1,0}. We emphasize here that 
the behavior of statefinder parameters in the figures is due to a selected choice of model 
parameters while it can change for other values of parameters. 
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Figure 1: The statefinder parameters for NADE are plotted for the parametric values 
= 1.1, n = 0.03, Cg = 0.02, C4 = 0.03, C5 = 0.04 and = i; = 1.1 




Figure 2: The statefinder parameters for HDE are plotted for the parametric values 
= 0.9, c = 0.5, Ci = 0.02, C2 = 0.3,and m = 1.1 and v = 1.2 
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